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ABSTRACT 

The reliability of astronomical observations at millimeter and sub-millimeter wave- 
lengths closely depends on a low vertical content of water vapour as well as on high 
atmospheric emission stability. Although Concordia station at Dome C (Antarctica) 
enjoys good observing conditions in this atmospheric spectral windows, as shown by 
preliminary site-testing campaigns at different bands and in, not always, time over- 
lapped periods, a dedicated instrument able to continuously determine atmospheric 
performance for a wide spectral range is not yet planned. In the absence of such mea- 
surements, in this paper we suggest a semi-empirical approach to perform an analysis 
of atmospheric transmission and emission at Dome C to compare the performance for 
7 photometric bands ranging from 100 GHz to 2 THz. Radiosoundings data provided 
by the Routine Meteorological Observations (RMO) Research Project at Concordia 
station are corrected by temperature and humidity errors and dry biases and then 
employed to feed ATM (Atmospheric Transmission at Microwaves) code to generate 
synthetic spectra in the wide spectral range from 100 GHz to 2 THz. This approach is 
attempted for the 2005-2007 dataset in order to check its feasibility. To quantify the 
atmospheric contribution in millimeter and sub- millimeter observations we are consid- 
ering several photometric bands, largely explored by ground based telescopes, in which 
atmospheric quantities are integrated. The observational capabilities of this site at all 
the selected spectral bands are analyzed considering monthly averaged transmissions 
joined to the corresponding fluctuations. Transmission and pwv statistics at Dome C 
derived by our semi-empirical approach are consistent with previous works. It is evi- 
dent the decreasing of the performance at high frequencies. We propose to introduce 
a new parameter to compare the quality of a site at different spectral bands, in terms 
of high transmission and emission stability, the Site Photometric Quality Ratio. The 
effect of the instrument filter bandwidth is involved on the estimate of the optical 
depth performed by the water vapour content knowledge. 

Key words: Site testing - Atmospheric effects - Submillimctre - Cosmology: obser- 
vations. 



1 INTRODUCTION 

Astronomy and astrophysics in the 100-1000 GHz band al- 
low the study of a large variety of processes, in the local and 
distant universe, which involve cool matter absorbing and re- 
radiating efficiently at these frequencies, in environments of- 
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ten unaccessed through observations at visible wavelengths. 
In fact, many key topics in modern astronomy and cosmol- 
ogy, such as galaxy formation and evolution, the amount and 
role of dark matter and dark energy in the universe, star 
formation, protoplanetary disks, or the properties of cold 
debris at the outskirts of the solar system, are related to ra- 
diative phenomena in this band. The field has undergone a 
huge development in the last two decades, thanks to the de- 
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Figure 1. Daily values of precipitable water vapour (pwv) estimated from the 12:00 UTC radiosounding measurements performed at 
Dome C over the period from May 2005 to January 2007 arc shown in the upper left panel. Monthly averages of pwv are overplotted 
as red dots. In the right panel the corresponding pwv vs. cumulative frequency is plotted. The bottom panel shows the monthly pwv 
fluctuations quantified as the daily values standard deviation, cr pwv . 



velopment of sensitive detectors, large cameras, polarization 
sensitive devices and spectroscopically capable instrumenta- 
tion. Some key achievements range from the measurement 
of the intensity and polarization power spectra of the cosmic 
microwave background at millimeter (mm) wavelengths, to 
the discovery and the characterization of the optically elu- 
sive sub-millimeter (sub-mm) galaxy population (SCUBA, 
BLAST), and the recent galaxy cluster surveys through sub- 
arcminute resolution observations of the Sunyaev-ZePdovich 
effect (ACT, SPT). 

Ground-based observations in the mm/sub-mm band are 
usually plagued by the transparency of the atmosphere (and 
its stability over time), mainly because of the presence of 
large, time-dependent pressure-broadened features in the 
emission (and absorption) spectrum of the water vapour. 
Of course, this issue is strongly mitigated when operat- 
ing stratospheric balloon-borne or airborne detectors, and 
completely averted when moving detectors on spacecrafts. 
BOOMERANG, BLAST, SOFIA, Planck and Herschel 
have proven the effectiveness of mm and sub-mm obser- 
vations from the stratosphere and from space, providing 
ground breaking advancements in their respective fields at 
the time of their operation. 

Anyway, the practical limitations on the telescope size, the 
weight and the accessibility of instrumentation still make 
substantially unfeasible the deployment of large (10m class) 
telescopes on balloons, aircrafts or satellites. As a matter of 
fact, the ground-based solution appears presently the only 
viable way to routinely perform high angular resolution ob- 
servations of compact objects and/or small spatial and spec- 
tral features in cool diffuse media at sub-mm wavelengths. 
As a consequence, the last few years have witnessed increas- 
ing efforts in the design and construction of large telescopes 
in places of the planet which provide the potentially most 



attractive atmospheric features for mm and sub-mm astron- 
omy. The community has realized the need to perform a 
thorough characterization of astronomical sites in terms of 
atmospheric opacity and stability across the whole mm/sub- 
mm spectral region, both for observation planning and for 
transparency monitoring during the observing sessions. 
At a time where bolometric detectors can easily approach 
the photon noise limit, and large cameras with hundreds 
or thousands of such detectors already allow to break this 
limitation, it is straightforward to realize that an improper 
characterization of the atmospheric properties may become 
the strongest restriction to the effective science return from 
ground-based instruments of the present (ACT, SPT) and 
next generation (CCAT). 

In order to continuously monitor the atmospheric transmis- 
sion several approaches are possible: tippers or tau-meters, 
hygrometers, GPS, water vapour radiometers, radiosound- 
ings and spectrometers. The first approaches allow a con- 
tinuous data recording by simple instruments but with the 
drawback of single frequency observations, needing a syn- 
thetic atmospheric model to infer transmission at other fre- 
quencies. 

Dome C is considered one of the best sites in the world 
to perform observations in a wide range of the electromag- 
netic spectrum allow ing a lso to explore T e raher tz windows 
jMinier et all (|2008h and iTremblin et all l|2()l \b). Anyway 
a wide frequency coverage transmission measurements cam- 
paign at Dome C, employing the direct spectroscopic infor- 
mation derived by an interferometric experiment, was never 
carried out. 

The goal of this paper is to compensate the lack of those data 
by estimating the atmospheric performance with a semi- 
empirical approach. We test this method using the avail- 
able dataset of radiosounding measurements recorded by 
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Figure 2. Daily values of precipitable water vapour (pwv) ob- 
tained through the present analysis from surface-level to 8 km 
amsl, and plotted versus W, the corresponding values of precip- 
itable water derived by Tomasi et al. (2011a) over the altitude 
range from surface-level to 12 km amsl. The data are best-fitted 
by a regression line with intercept equal to 9.46610 - 3 and slope 
coefficient equal to 0.9425, which was obtained with regression co- 
efficient R = +0.993, and provided a standard error of estimate 
SEE = 0.04 mm. 
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Figure 3. Atmospheric transmission spectra as modeled by ATM 
program for each radiosounding. Photometric bands in Table [T] 
(gray) match the main transmission windows. 



Table 1. Characteristic spectral bands assumed in this work. 





v (GHz) 


A (Aim) 


FWHM(%) 


References 




150 


2000 


22 


1, 2, 3, 4 


LF 


220 


1400 


13 


1, 2, 3 




270 


1100 


18 


1, 2, 3 




350 


860 


8 


1, 2, 5, 6 




660 


450 


11 


5, 6 


HF 


870 


350 


13 


5 




1500 


200 


5 


7 



the Routine Meteorological Observations (RMO) Research 
Project (www.climantartide.it) at Concordia station in the 
period from May 2005 to January 2007, carefully corrected 
for the main lag errors and dry biases. The profiles of air 
temperature, pressure and relative humidity allow to gener- 
ate synthetic sp ectra, ranging from 100 GHz to 2 THz, with 
the ATM code jPardo et all (|2001al )). 

The paper is organized as follows. Atmospheric synthetic 
spectra, as derived by ATM code, are described in Sect. [2] 
In Sect. [3] estimates of atmospheric transmission and emis- 
sion corresponding to largely explored ground based tele- 
scope bands between 150 and 1500 GHz are analyzed. The 
effect of the filter bandwidths on the estimate of opacity is 
for the first time included in the relation showing a contri- 
bution up to a 30 per cent over-estimate on the opacity in 
the case of the highest frequency band. 
A parameter to rank the observational conditions for each 
of the selected spectral bands is introduced as the ratio be- 
tween average transmission and the corresponding fluctua- 
tions. 

Finally a discussion on the analysis and the conclusions are 
summarized in Sect. [4] 

A detailed description of the correction procedure used to 
analyse the raw radiosounding data and determine the ver- 
tical profiles of the main thermodynamic parameters is re- 
ported in the Appendix. 



References: (1 ) SPT , ISchaffer et "ahl J2011I1; (2} MITO, 

iDe Petris et al.l [120021) ; ( 3) A CT, ISwetz et alJ d201lll: (4) 
BRAI N, iBattistelli etafl J2012h; (5) SCU BA, Irlolland et alj 
dl999h: (&) SCUBA-2 iDemnsev et aT] l|201fj l: (7) THUMPER, 
IWard-Thompson et afl d2005h ; 

2 SYNTHETIC SPECTRA PRODUCTION 

At present, for the site of Dome C we can rely only on the 
atmospheric monitoring performed at a few individual fre- 
quencies, with no simultaneous measurements in different 
regions of the spectrum. In order to compensate for the lack 
of a continuous and spectrally wide atmospheric monitoring 
at Dome C, we predict the performance in the mm/sub- 
mm bands in the period from May 2005 to January 2007 by 
means of a semi-empirical approach. 

A set of raw radiosounding data was recorded for the present 
study, consisting of an overall number of 469 radiosounding 
measurements taken routinely at Dome C, at 12:00 UTC 
from May 2, 2005 to January 31, 2007 ranging from a min- 
imum of 15 in May 2005 to a maximum of 30 in November 
2006. 

In general, each radiosonde measurement consists of values 
of air pressure P, air temperature T and relative humidity 
RH, taken at more than 800 standard and additional levels 
in the altitude range from surface to 10 km above mean sea 
level (amsl). Data provided by the radiosonde sensors are 
affected by lag and instrumental errors as well as by various 
dry biases. They were all corrected following the procedure 
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described in the Appendix. 

The time-patterns of the daily pwv values are shown in Fig. 
[1] Two main features are evident in Fig. [T] showing that 
the majority of pwv values are lower than 0.3 mm during 
the austral autumn months, although presenting largely dis- 
persed patterns (of 50 per cent or more), and, hence, low 
stability. 

As shown by iTomasi et al.l l|2011al ). a limited contribution 
is given to the overall value of atmospheric pwv by the 
amount of water vapour present in the Upper Troposphere 
and Low Stratosphere (UTLS) region from 8 to 12 km amsl, 
while negligible fractions of pwv ranging mainly between 
0.003 and 0.005 mm throughout the year are present in the 
stratosp here from 12 to 50 km above Dome C (jTomasi et al.l 
ll2011blk 

To verify the reliability of the present estimates of pwv, 
a comparison is made in Fig. [2] among the present daily 
values of pw v and those correspondingly determined by 
ITomasi et al.l (|2011al ) (indicated as W) using a more ad- 
vanced correction procedure from surface-level to 12 km 
amsl. The comparison showed that a close re lationship ex- 
ists bet ween the present results and those of ITomasi et alj 
|2011al 'l. defined by a regression line with nearly null inter- 
cept and slope coefficient of + 0.9425, having regression co- 
efficient better than + 0.99, and providing a standard error 
of estimate equal to 0.04 mm. These findings clearly indi- 
cate that the present evaluations of pwv, as obtained over 
the altitude range from surface-level to 8 km amsl, are fully 
suitable for the purposes of our study, especially considering 
the intrinsic uncertainty of the simulation model. 

We have estimated synthetic spectra in emission 
and in opacity by means of the ATM code in the wide 
spectral range from 100 GHz to 2 THz. Each spectrum is 
derived considering the corrected radiosounding data. The 
transmission corresponding to each radiosounding dataset, 
estimated from optical depth spectra as T = e~ T , is shown 
in Fig. H 

In the period under consideration, the inferred pwv values 
show an average close to 0.3 mm with a mean dispersion 
of about 150 (im (see Fig. [IJ. The same amount of pwv 
variation can contribute with a different weight to the total 
optical depth. As example in Fig. [4] we represent the optical 
depth fluctuations derived by ATM, quantified as the max- 
imum dispersion, due to fluctuations of pwv of the order 
of 150 (jm around three different pwv values (0.15, 0.5 and 
1.0 mm). It is worthy of note that for low pwv content, At 
can be as high as 60 per cent in the high frequency windows. 



3 MULTI-BAND ANALYSIS 

A quantitative analysis is performed considering 7 photo- 
metric bands centered at the frequencies of several astro- 
physical and cosmological experiments: South Pole Tele- 
scope (SPT), Atacama Cosmology Telescope (ACT), Mil- 
limetre and Infrared Testagrigia Observatory (MITO) and 
BRAIN (B-mode RAdiation INterferometer) for Low Fre- 
quency (LF) atmospheric windows; Submillimetre Common- 
User Bolometer Array (SCUBA and SCUBA-2) and Two 
HUndred Micron PhotometER (THUMPER) for sub-mm 
bands (High Frequency, HF). The central frequency of each 
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Figure 4. Optical depth fluctuations corresponding to a 150 fim 
variation around three selected pwv values: 1 mm (red line), 500 
jim (green line), 150 fim (black line). 

band, as well as the bandwidth, quantified with the FWHM 
(Full Width Half Maximum) , are listed in Table Q] (see also 
Fig. [3} . The band profiles are assumed to be top-hat assum- 
ing in this way the maximum rejection to off-band contri- 
butions. 

To assess the constraints on astronomical observations 
arising from the atmosphere emission above Dome C, we 
give an estimate of the NEP (Noise Equivalent Power) and 
the NEFD (Noise Equivalent Flux Density) for all the seven 
bands. In fact in such a wide spectral region both the quan- 
tities are normally employed: the power density, mainly for 
the low frequency bands, while the flux density, for the high 
frequency region. The quoted NEP is the root of the sum 
of NEPatm, the term considering the atmospheric emission 
fluctuations, and NEP^ ele , i.e. the instrumental contribu- 
tion to the photon noise. The atmospheric emissivity spec- 
tra are generated by ATM. The telescope is assumed a 10-m 
in diameter Al-mirror with a surface emissivity of the or- 
der 3 per cent at 150 GHz and depending on the frequency 
as yjv. The throughput of the telescope is assumed diffrac- 
tion limited at each band. Focal plane optical efficiencies are 
taken as unitary for all the bands as well as telescope main 
beam efficiency. The dominant sky sources (CMB and dust) 
are not included, the instrument detector noise is assumed 
lower than the background noise and the spillover emission 
is neglected. In order to quantify the maximum variation of 
these quantities we plot in Fig. \E\ NEP and NEFD val- 
ues for all the bands, for the extreme conditions occurred 
during the austral summers and winters at Dome C in the 
2005-2007 period. 

3.1 Dome C statistics comparison 

To validate the proposed semi-empirical approach, we 
compare the derived atmospheric performance with the 
results available in literature. 

Fig. [6] shows Dome C atmospheric transmission as a 
function of the cumulative time frequency derived by 
radiosounding data and ATM model for the bands listed 
in Table [1] (the corresponding quartiles are reported 
in Table |2). Transmission statistics at Dome C per- 
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Figure 5. Noise Equivalent Power (upper panel) for the seven 
bands in two extreme atmospheric conditions in the austral winter 
(red triangles) and summer (black diamonds). In the lower panel, 
the same for the Noise Equivalent Flux Density. 



Table 2. Transmission quartiles matching cumulative distribu- 
tions in Fig. [6] 



vo (GHz) 


A (Aim) 


25% 


50% 


75% 


150 


2000 


0.97 


0.97 


0.96 


LF 220 


1400 


0.95 


0.94 


0.94 


270 


1100 


0.94 


0.93 


0.92 


350 


860 


0.91 


0.89 


0.87 


660 


450 


0.64 


0.56 


0.46 


HF 870 


350 


0.58 


0.50 


0.40 


1500 


200 


0.15 


0.08 


0.03 


Table 3. Pwv quartiles comparison. 


Period 


25% 


50% 


75% 


-teferences 


01/1997 


0.38 


0.52 


0.68 


1 


05/2005-01/2007 


0.20 


0.30 


0.45 


2 


2008 


0.15 


0.24 




3 


2008-2010 


0.21 


0.27 


0.35 


4 


12/2009-01/2010 


0.49 


0.75 


1.1 


5 



References: (11 IValenziano fc DaU'Oglid fll999T) ; (2) this work; (3) 
lYang et al] (|201Cll l: C4) iTremblin et al.l l|201 J i~f5'l lBattistelli et all 
||2012|1 . 



1.00 



0.95 



E 0.90 



0.85 



0.80 
0.80 



0.60 



0.40 



0.20 



LF Bands 




0.00 



350 GHz (0.86 mm) 
270 GHz (1.1 mm) 

220 GHz (1.4 mm; 
150 GHz (2.0 mm; 



HF Bands 
(1500 GHz) 




0.0 



0.2 



0.4 0.6 
Frequency T 



0.8 



Figure 6. Atmospheric transmission vs. cumulative time fre- 
quency for Dome C corresponding to the atmospheric windows 
listed in Table [T] 



forme d by IValenziano fc Dall'Oglid ill999l). iMinier et all 



J2008h. lYanget alJ (j2010T l. ITremblin et all (|201ll ) and 



iBattistelli et all (|2012h are compared with our analysis. 
Low frequency atmospheric windows show high trans- 
parency during the whole period confirming that high 
quality mm observations can be performed from this site 
for most of the time. For instance the 150 GHz 50 per cent 
quartile transmission is about 97 per cent (see the green 
line in Fig. [6jl . This is consistent with the 95 p er cent value 
recently measured by IBattistelli et~ai1 \2QV& \ during the 
summer campaign in December 2009/ January 2010, even 
considering their integrated in-band result. 
Median transparency for the 220 GHz atmospheric window 
is about 95 per cent (see the cyan line in Fi g. [6]) as 
already derived bv IValenziano fc DairOgliol (|l999h by pwv 
measurements performed with a portable photometer in 
January 1997. 

Dome C 450 /im window remains above a tr ansmission of 
60 pe r cent for 50 per cent of the time. In IMinier et al.1 
(2008) the atmospheric transmission at Dome C has been 
calculated throug h the 5-years pwv da ta from the South 
Pole available in IPeterson et al.l (|2003l ) and extrapolating 
the corresponding at mospheric transm ission at Dome C 
using the model in lLawrencd l|2004h . They found that 
450 [im median transmissio n at D ome C is about 70 per 
cent. Recently lYang et al.1 l|2010l) measured a 450 fim 
median winter transmission at Dome C of about 60 per 
cent estimating pwv with the Microwave Humidity Sounder 
(MHS) sounding on the National Oceanic and Atmospheric 
Administration (NOAA) ozonesondeas in 2008. 
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Figure 7. Monthly average values of transmission (T) and the relative monthly fluctuations or, plotted as rms values, shown in different 
colors: austral winter months in the upper panel and summer months in the bottom panel. Monthly Site Photometric Quality Ratio 
(SPQR) values are shown in the right panel. 



Dome C median transmission for the 350 \im atmo- 
spheric window is about 5 per cent (see t he bo ttom panel 
in Fig. [6]), as derived by iTremblin et all (|201lf ) using the 
MOLIERE model and 200 fim optical depth measurements. 
They found also that the Dome C 200 fim window opens 
with a transmissi on of 10 per cent fo r less then 25 per cent 
of the time while lYang et al.l (|2010T l found that the trans- 
mission at 200 fim is about 13 per cent for 25 per cent of 
the time in 2008. 

The 200 (jm transmission as a function of the cumulative fre- 
quency is the black solid line in the bottom panel of Fig(6] 
the 25 per cent quartile transmission value is above 10 per 
cent. 

Pwv quartiles since May 2005 until January 2007 (see the 
right panel in Fig. [TJ) have been compared with Dome C wa- 
ter vapour estimates performed in previous works in Table 

El 



3.2 High transmission and emission stability 

Following an observational approach, we report the statis- 
tics of integrated in-band quantities, like emission and 
transmission. Both monthly averages (T) and relative 
dispersions ctt (rms values) of in-band transmissions, 
ranging from May 2005 until January 2007 and splitting 
between austral summer (from October to February) and 
winter months (from March to September) are shown in 

Fig.m 

Monthly averaged transmission fluctuation is a good 
proxy of emission stability due to the fact that transmission 
and emission fluctuations are linearly correlated. In addi- 
tion we assume that the estimated monthly averaged fluctu- 
ations, quantified in terms of the standard deviation, could 



Table 4. Seasonal averages of the SPQR. 





u (GHz) 


Summer 


Winter 




150 


272 


335 


LF 


220 


127 


152 




270 


79 


94 




350 


36 


42 




660 


6 


7 


HF 


870 


5 


6 




1500 


1 


2 



be an underestimate of atmospheric stability because they 
derive from a daily data sampling, the time interval between 
two consecutive radiosoundings. 

We note that during the austral winter the atmospheric 
transmission in all the considered bands is generally higher, 
as expected. (T) shows values close to the unity in mm bands 
and decreases towards THz windows, while relative disper- 
sions <tt have the opposite trend. As an example, the best 
atmospheric conditions (in term of high transmission) occur 
when the atmospheric fluctuations ot are larger than oth- 
ers months (red dots in Fig. [7}. Referring to the atmospheric 
window centered at 200 /Ltm, when the transmission has the 
maximum value, the large fluctuations at short time scales 
are likely to degrade the quality of a scientific observation. 
In addition it is not possible to identify the month with 
the best atmospheric performance as one can see from the 
gap between two consecutive years atmospheric transmis- 
sion and fluctuations (red and blue dots in Fig. [TJl . 
450/xm and 350/im bands transmission show a reduction of 
few percent ranging from winter to summer months, while 
fluctuations are not sensitive to seasonal effects. 
All the considered bands are characterized by high stability 
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in October (see the black or orange dots in the middle panel 
of Fig. [7)l with the exception of the 200 /im window, showing 
high stability especially during summer months like January 
or February, when atmospheric transparency is not suitable 
to perform astrophysical observations. 

To quantify the real capability of the observational site we 
need to study the atmospheric performance, mainly the 
stability, strongly affected by the weak reproducibility of 
weather conditions at long time scales. In order to highlight 
this issue we introduce a specific parameter, the Site Pho- 
tometric Quality Ratio: 

SPQR = @1 (1) 

relating monthly averaged transmission to its fluctua- 
tions, sampled on a daily timescale, for all the considered 
atmospheric windows. SPQR amplitude provides informa- 
tion about atmospheric performance and it allows us to der- 
mine if high transmission combined with high transmission 
(i.e. emission) stability conditions are both satisfied for each 
band. Even if we are not able to identify the desired SPQR 
threshold, this factor could represent a useful tool to com- 
pare several bands performance or sites. In the right panel 
of Fig.[7]monthly values of the Site Photometric Quality Ra- 
tio are shown in different colors. The differences between the 
two years are more evident in SPQR, anyway a decrement of 
the Site Photometric Quality Ratio towards THz regime oc- 
curs in austral winter as well as in summer periods. Seasonal 
averaged values of the Site Photometric Quality Ratio in Ta- 
ble 2] suggest the good quality of atmospheric conditions in 
the low frequency bands, notably during the austral winter. 
While SPQR appears useful for comparison among different 
bands at Dome C it could also be employed for comparison 
among different sites. It is worth reminding that it is diffi- 
cult to quantify for SPQR a threshold value to discriminate 
the goodness of a site. 

Two outcomes can be gathered from this analysis. If we be- 
lieve in the transmission values, as derived by this semi- 
empirical approach, a continuous atmospheric sampling is 
mandatory at least at high frequency to contrast the low 
transmission stability. Otherwise if we consider the data 
derived in this work not reliable enough for an accurate 
estimate of atmospheric properties, we need direct obser- 
vational techniques. In either case continuous atmospheric 
transparency measurements in all the spectral range of in- 
terest are necessary. 



3.3 Effect of broadband filter on optical depth 
estimate 

The average of the optical depth over a band, t„ (Az/), is 
larger than its central value t vq so the opacity is overesti- 
mated by broa dband instruments l ike tippers, as remarked 
as example by ICalisse et al.l (|2004T l . The determination of 
this effect is not unique because several pwv values could 
give the same in-band integrated opacity. Low-frequency in- 
struments are less sensitive to this degeneracy even for large 
values of the bandwidth due to the flatness of the corre- 
sponding atmospheric windows. On the other hand a sub- 
mm broadband instrument overestimates the opacity (un- 
derestimates the transmission) and this difference depends 




0.2 0.4 0.6 0.8 1.0 1.2 

PWV (mm) 



Figure 8. Best fit of the correlation between the narrow-band 
opacity data and atmospheric water vapour for mm windows (top 
panel) and for the sub-mm bands (bottom panel). 



on the filter shape as well as on the the atmospheric con- 
ditions. Little variations of atmospheric conditions give rise 
to a dispersion of this overestimate because of the relative 
shapes of the atmospheric window and the corresponding fil- 
ter. For each band in Table [T] we have included the effect in 
the relation between the integrated zenith opacity t vq (Au) 
and pwv values the effect of the instrumental bandwidth Av 
(see Fig. [8] and Table 0: 

t V0 (Av) = (o + a-^Av) + (b + ^A^pwv (2) 

ao and bo are the linear fit coefficients of the t„o vs 
pwv relation referred to a narrow band filter matched to 
the central frequency and ai and bi take into account the 
dependency on the instrumental bandwidth Av, linearly ap- 
proximated at least in the range within the maximum band- 
widths as reported in Table [T] Realistic band profiles could 
highlight the effect instead of our approximation with top- 
hat profiles. The net result is that the optical depth can be 
overestimated at most of 30 per cent at 200 \im, assuming 
the pwv best quartile from Table[3] while low frequency win- 
dows are less sensitive to this effect, as expected (10 per cent 
at 150 GHz). 

The uncertainty related to the optical depth value due to the 
intrinsic scatter of the to vs pwv relation, can be approxi- 
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Table 5. Opacity-pi™ relation best fit parameters evaluated for the bands of interest. 





vo (GHz) 


ao 


Ol 


bo (mm 1 ) 


6i 


CO 


Cl 




150 


0.012 


0.0077 


0.024 


0.0054 


0.0016 


0.0012 


LF 


220 


0.015 


0.052 


0.053 


0.030 


0.0024 


0.0012 




270 


0.024 


0.016 


0.088 


0.0070 


0.0040 


0.0016 




360 


0.032 


0.073 


0.19 


0.41 


0.0093 


0.014 




660 


0.16 


0.021 


1.35 


1.71 


0.066 


0.073 


HF 


870 


0.18 


0.68 


1.51 


1.96 


0.074 


0.070 




1500 


0.53 


10.13 


5.82 


12.69 


0.29 


0.35 



mated by a linear trend as a function of the instrumental 
bandwidth: 

a Tuo (Ai/) = co + ci Av (3) 

The optical depth uncertainty turns out to be 0.002 at 
150 GHz and rise up to 0.3 at 200 fj,m, assuming the dis- 
persion independent on pwv value (see Fig. [SJ. As a conse- 
quence the percentage uncertainty on optical depth estimate 
is about 15 per cent all over the considered atmospheric win- 
dows assuming the best pwv quartile and it remains above 
10 per cent even assuming the 75 per cent quartile in Table 

El 

The six parameters corresponding to the seven bands are 
listed in Table [5] The Eq. [2] is useful to infer the atmo- 
spheric opacity at the preferred frequency, with a specific 
bandwidth, when the pwv content is known, but it is impor- 
tant to remind that this relation is appropriate only in the 
en virons of Dome C. 

In lTremblin et all (|20 1 if ) the opacity is related to the atmo- 
spheric pwv by means of the MOLIERE model. The result- 
ing linear regression of the pwv as a function of the 200 fim 
opacity and the corresponding best fit parameters in Table 
[5] neglecting a± and 6i, gives less than 5 per cent difference 
in transmission for low pwv values. Such a gap could be eas- 
ily included in the atmospheric performance variations ob- 
served at Dome C over the years. Also the difference in trans- 
mission evaluated for 220 GHz bes t fit parameters in Table 
!5iand to ( 22HGH z )-pwv linear fit in IValenziano fc Dall'Oglicl 
( 1999) is lower then 4 per cent. 



4 CONCLUSIONS 

The quality of cosmological and astrophysical measurements 
performed from ground based observational sites in the mm 
and sub-mm wavelength regions are strongly dictated by 
the atmospheric performance. 

The simultaneous measurement of atmospheric trans- 
parency and transmission fluctuations, i.e. emission 
stability, is a necessary condition to determine the true 
capabilities of the site of interest. 

We try to monitor the atmosphere across a wide spectral 
range, mm and sub-mm, with a semi-empirical approach. 
The transmission at Dome C is inferred by generating 
ATM synthetic spectra as derived by radiosounding data 
in the period from May 2005 to January 2007. Excellent 
performance is evident in the low frequency bands while 
large emission fluctuations are present in the high frequency 
bands. In fact even if the median winter transmission is 



large in all the considered atmospheric windows, daily 
atmospheric emission fluctuations are not negligible and 
become remarkable in the sub-mm range. In addition, large 
timescales fluctuations of the atmospheric performance 
have been detected during two consecutive years. 
The ratio between monthly averaged transmission and 
the corresponding fluctuations, defined Site Photometric 
Quality Ratio, turns out to be an efficient estimator to rank 
the photometric performance of the atmosphere, in terms 
of stability, above Dome C, as well as any observational 
site. It allows to verify when high transmission as well as 
low skynoise requirements are satisfied for the atmospheric 
window of interest. The SPQR threshold for each band 
is not easily defined: it is depending on the detectors 
architecture and on the adopted observational strategy. 
We attempted to validate the proposed semi-empirical 
approach comparing pwv and transmission quartiles with 
other site-testing campaigns performed at Dome C during 
the last years also at different wavelengths. 
In the usual linearly dependent opacity-pu™ relation, we 
include the effect due to the bandwidth of the monitor 
instrument. 

Anyway only direct and frequent measurements of atmo- 
spheric transmission in a wide spectral range can provide a 
perfect knowledge of atmospheric influence on astronomical 
observations. If the opacity measurements are done in a 
narrow (a few MHz) spectral coverage, it is impossible 
to distinguish between clear sky opacity, hydrometeors 
contributions, and systematic errors. A wide frequency 
coverage (several hundreds of GHz) is necessary to make 
sure we are in clear sky conditions and no instrumental 
offset is affecting our measurement and our analysis. In this 
way it is also possib le to determ i ne the dry and the wet 
continuum terms, see iPardo et al.l l|2001bl ). 
A large spectral sampling can be achieved at the price 
of a bit complex instrument. The possibility to monitor 
the atmosphere towards different positions in the sky, also 
avoids bias due to a spatial model assuming the multi layers 
approximation. 

A dedicated spectromete r, like the one proposed for Dome 
C (|De Petris et al.l ([20051) 1 and in operation at Testa Grigia 
station (3500 m a.s.L, Alps, Italy) in a spectrally li mited 
version (100 -5- 360 GHz), CASPER 2 (jPecina et al.l (|2010h 
and De Petris et al. in prep.), is a viable solution. 
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APPENDIX A: CORRECTION OF THE 
RADIOSONDE DATA AND CALCULATIONS 
OF PRECIPITABLE WATER VAPOUR. 

The meteorological data were obtained at Dome C using 
two Vaisala radiosonde models: (i) the RS92 model for 430 
measurement days, i.e. for 94 per cent of the overall days, 
and (ii) the RS80-A model for 29 radiosonde launches only. 
Each triplet of signals giving the measurements of P, T and 
RH at a certain level was sent by the transmitter onboard 
the radiosonde to the ground station every 2 s. Considering 
that the radiosonde ascent rate was in general 5 - 6 m s _1 , 
the triplets of signals were recorded in altitude steps of 10 - 
12 m. 

The main characteristics of the three sensors (Barocap, 
Thermocap, and Humicap) mounted on the two rad i osond e 
models are available in Table 1 of iTomasi et al. I i|2006l ). 
where their measurement range, resolution, accuracy, re- 
peatability in calibration, and reproducibility in sounding 
are given. 

The measurements of P, T and RH provided by the ra- 
diosonde s ensors were al l corre cted following the procedure 
defined bv lTomasi et al.l l|2006l ). which consists of numerous 
steps adopted to minimize the errors due to: 

(i) the not correct calibration of the Barocap sensors; 

(ii) the effects caused by solar and infrared radiation heat- 
ing, heat conduction and ventilation on the Thermocap sen- 
sors; 

(iii) lag errors, ground-check errors, and dry biases of the 
Humicap sensors due to basic calibration model, chemical 
contamination, tempe rature dependence and sensor aging, 
corrected according to lWang et al.l (2002). 

This procedure subs tantially differs from that defined by 
ITomasi et alj (|2011al ) only in the parts regarding the correc- 
tion of solar heating dry biases for both A- and H-Humicap 
se nsors: (1) those of th e A-Humicap sensor were corrected 
by Tomasi et al. 1 (l2006l) usi ng the algorithm of iTurner et al.l 
(2003). while lTomasi et all ()2011al) preferred to use the algo - 
rithm derived more recently by Cadv-Pereira et all |2008); 
and (2) those of the H-Humicap sensor were corrected by 



ITomasi et al. (l2006h using the avera ge correction factors 
proposed by Miloshevich et all (l2006h as a function of so- 
lar zenith angle, while ITomasi et al" l|2011al ) employed the 



pair of day-time and n ight-time correction algorithms of 
iMiloshevich et all |2009). A large part of the few percent 
discrepancies found in the comparison shown in Fig. [3] be- 
tween the present va lues of pwv and those determined by 
ITomasi et ail (|2011al ) arise from the use of these different 
algorithms in correcting the instrumental and solar heating 
dry biases affecting the field measurements of RH. 

Using the correction procedures previously described, 
the daily vertical profiles of pressure P(z), temperature 



T(z) and RH(z) were determined at fixed levels above the 
surface-level, in regular steps of 25 m from 3.25 to 4 km, 50 
m from 4 to 5 km, and 100 m from 5 to 8 km amsl. 
In order to calculate the values of absolute humidity q(z) at 
the same fixed levels, the following procedure was adopted, 
consisting of the three steps: (i) calculation at each level of 
the saturation vapour pressure E(T) in the pure ph ase over 
a plan e surface of pure water, using the well-known iBoltonl 
( 1980) formula; (ii) calculation at each level of the water 
vapour partial pressure e(z) as the product E(T) x RH(z); 
(iii) calculation at each level of absolute humidity q(z) (mea- 
sured in g m -3 ) in terms of the well-known equation of state 
of water vapour, and, hence, as the ratio between e(z) (in 
hPa) and the product R w x T(z) (in K), in which the water 
vapour gas constant R w = 0.4615 J g~ l K~ x is put in place 
of the constant R a = 0.287 J g~ x K" 1 used in the equation 
of state for dry air. 

For all the 469 daily vertical profiles of q{z) obtained 
using the above procedure, the values of pwv were then cal- 
culated by integrating each vertical profile of q(z) from the 
surface-level to 8 km amsl (i.e. up to 4.767 km above the 
ground level). 
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